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' Abstract. We study vortex-like solutions to the Lifshitz-Chern-Simons theory. We 

' nn d that such solutions exists and have a logarithmically divergent energy, which 

suggests that a Kostclitz-Thouless transition may occur, in which voxtex-antivortex 
pairs are created above a critical temperature. Following a suggestion made by Callan 
and Wilzcek for the global U(l) scalar field model, we study vortex solutions of the 
Lifshitz-Chern-Simons model formulated on the hyperbolic plane, finding that, as 
expected, the resulting configurations have finite energy. For completeness, we also 



o 



explore Lifshitz-Chern-Simons vortex solutions on the sphere. 
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1. Introduction 

First proposed in [TJ, Lifshitz-Chern-Simons (LCS) theory can be understood as a 
modification of the 2+1 dimensional z = 2 Lifshitz scalar theory [2], [3] by the addition 
of a non-local term. The action is obtained by dualizing the scalar field into a gauge 
vector and then deforming the resulting Lagrangian with a Chern-Simons three-form. 
In terms of the original scalar field, the Chern-Simons term corresponds to a non-local 
deformation. 

As originally presented, the Lifshitz-Chern-Simons theory models a system that 
experiences an isotropic to anisotropic phase transition at zero temperature. In the 
anisotropic phase, the £0(2) rotational symmetry enjoyed by the action is spontaneously 
broken by the electric field acquiring a non-zero vacuum expectation value. 

In [3], the LCS action was shown to be equivalent to the model introduced in 
[5], [6] to describe the low energy behavior of a charged spinless 2d fluid in the presence 
of an external perpendicular magnetic field. The anisotropic phase appears through 
renormalization of the free parameters of the theory, and reproduces the phenomenology 
in the presence of a parallel magnetic field [7] . 

In the present work, we study vortex solutions of the Lifshitz-Chern-Simons 
action. We find that such solutions exist and have a logarithmically divergent energy. 
An heuristic argument then suggest that they may be entropically favored at high 
temperatures, leading to a Kosterlitz-Thoules like transition. 

It was suggested in [8] that the logarithmically divergent energy of the global vortex 
solution to the charged scalar field model, could be cured by formulating the model in a 
negatively curved space. Following such suggestion we study vortex solutions to the LCS 
model in the hyperbolic plane. For completeness we also analyze the vortex solution on 
the sphere. We find that the resulting configurations have finite energy. At this point, 
one may wonder whether it makes sense to formulate the theory in curved space. We 
would like to point out that a scalar theory in curved space, that reduces to a Lifshitz 
scalar with z = 2 on a flat metric, turns out to describe the off-plane fluctuations of a 
tensionless membrane [15j . 

The plan of the paper is the following: in section II we present the Lifshitz-Chern- 
Simons theory defined on a general curved manifold. In section III we show that vortex 
solutions exist in this model by solving numerically the Euler-Lagrange equations, the 
relevant properties of the solution are also discussed in this section. In section IV we 
give a brief summary of our results. In the appendix A we review the relation between 
the Lifshitz-Chern-Simons theory and the Lifshitz scalar one. 
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2. The Lifshitz-Chern-Simons action 

The Lifshitz-Chern-Simons action in an arbitrary 2-dimensional curved space with 
metric is defined as 

S = J dtd 2 x^g- [e^dta, - d t a t ) - ^ (/^(V^) 2 + 6 2 ) - 

+ k \aj> - ^4 a ^ a iJ " \ (m 2 ( e i) 2 + ^( e *) 4 )) > W 

where repeated (or squared) lower latin indexes are understood to be contracted with 
the curved inverse metric g % \ i.e. Xii/i = g^Xii/j. V, is the standard covariant derivative, 
and b is a shorthand for b = e 1 -* diCij j \fg with e 4 - 7 the Levi-Civita symbol e 00 = e 11 = 
and e 01 = — e 10 = 1. The couplings k 2 , k, A are dimensionless, while [m 2 ] = 1/L 2 . Stated 
differently, the Weyl (scale) transformation of the metric — > Q 2 gij is a symmetry of 
the first three parentheses and the A-deformation provided one simultaneously scales 
the time coordinates as t — > Q 2 t, a t — > Q~ 2 a t . 

By electric/magnetic duality the first two parentheses are mapped to the Lifshitz 
scalar action, while the third one, representing a Chern-Simons deformation, leads to a 



non-local term for the Lifshitz scalar field (see appendix Appendix A ). Finally, the last 
parentheses in (JT|) takes into account possible deformations of the theory. 
The equations of motion derived from the above action read 

Vi e 4 + kb = , (2) 
d t ai - didt + K 2 V 2 ei - (m 2 + A(ej) 2 ) e, L = , (3) 

g ld d t e j + — (djb + k(d t aj - dja t )) = , (4) 
v 9 

where V 2 = g^ViVj. The first equation can be used to obtain the magnetic field b once 
the electric field et is known. The second equation on the other hand defines the electric 
field in terms of the vector potentials a t ,cij. Combining the first and second equations 
to eliminate b and at, aj from the third, we find 

kg^dtej = ^- (VjVfce* + k 2 (Vv 2 ^- - (m 2 + A(e fc ) 2 ) e,)) . (5) 

Since we are interested in static solutions, we will drop the time derivative appearing 
on the left hand side. In the following sections we will solve for using ([5]) and obtain 
b from (ED. 

For static configurations, the energy functional takes the form 

E = \J d 2 x^-g ^ 2 (V^) 2 + ^(V,e 4 ) 2 + m 2 (e,) 2 + ^(e,) 4 + e\ , (6) 

where we have added a zero point (constant) contribution E Q , which will be adjusted 
below. 
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3. Solutions 

3.1. Flat space 

In flat space = 5ij and the Weyl rescalling mentioned in the previous section can 
be realized by the coordinate transformation x l — > Qx l : the combined transformation 
(t, x l ) — > (Q 2 t,Qx l ) and (a^a^ej) — > (Q~ 2 a t , fi -1 ^, f£ _1 ej) is a symmetry of the first 
three parentheses of action (OQ). We therefore conclude that the first three parentheses 
in (CQ) describe a critical point with z = 2 dynamical scaling exponent. 

The last parentheses in ([!]) amounts to possible local deformations of the Lifshitz 



action (see appendix Appendix A to see its expressions in terms of the Lifshitz scalar 
0). It is immediate to see that the term m 2 (e^) 2 is relevant, dominating in the IR. When 
m 2 > the {diej) 2 can be ignored and the theory flows to a z — 1 infrared fixed point. 
If m 2 < 0, a (e 2 ) 2 term with positive coefficient A becomes mandatory to stabilize the 
theory. In this situation, a classical expectation value for will appear (see [16], [17] for 
related computations involving the A-term). 

As discussed above, the order parameter of the theory is a vector on the 2- 
dimensional plane. We shall classify the topological disjoint classes of solutions by the 
winding of the vector on the circle at infinity. Explicitly, for solutions approaching 
infinity as (ei,e2) — > (cos(n9), sin(n#)) eo, the topological charge can be defined as 

n = — I e y e; dej , (7) 

where C is the circle at infinity and e, is the electric field normalized by its vacuum 
expectation value. Note that n takes integer values. 

3.1.1. Vacuum solution: The true vacuum of the theory corresponds to the lowest 
energy static solution of fl3])-dl]) and its symmetry will depend on the sign of m 2 . For 
m 2 > the minimum energy solution to (|5]) is = 0, we call this the "isotropic" phase. 
For m 2 < 0, a uniform vacuum expectation value develops = \j—m 2 j\ Ui, with Ui 
an arbitrary unit constant vector. This solution breaks the global SO (2) symmetry 
enjoyed by the flat space action (0Q) and we name it the "anisotropic" phase. It has 
n = winding number, and vanishing energy provided E = Am 4 /2. 

3.1.2. Vortex solution: In what follows we will consider solutions with n = 1 for the 
m 2 = — |m| 2 < case and write the flat metric in polar coordinates as 

ds 2 = dr 2 + r 2 d9 2 . (8) 

To simplify the equations it is useful to rescale r = R D p where R Q = a/1 + k 2 ^ 2 R a with 
R a = l/\m\k. The n = 1 static vortex ansatz corresponds to 

e P (p,0) = ^ f{p) , ee(p,9) = 0. (9) 
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Plugging it into <^ the resulting equation of motion reads 

f" + -f'-\f-f(f 2 -i) = 0- (10) 

Remarkably, this equation coincides with the relativistic n — 1 global U(l) vortex 
equation [10J. The appropriate boundary conditions for a vortex configuration are 

/(p)->0, p^O, (11) 
/(p)^l, p^oo. (12) 

A power series expansion close to the origin shows two independent behaviors f(p) ~ 
p ±x , the linear one being the proper choice for a non-singular vortex, 

f(p) = (3p + 0(p 3 ), p«l. (13) 

At infinity the radial profile asymptotes its vacuum expectation value as 



2p 2 



P > 1, 



(14) 



in this last expression we have dropped the exponentially decaying homogeneous piece 
that arises upon linearizing ( II Op at infinity. 

We numerically explored the space of solutions of eq. ffTO]) shooting from the origin 
for different values of (3 looking for the asymptotics ( IT41) . The solution we found is 
plotted in fig. (JTJ) with the slope of / at the origin being j3 ~ 0.58319. 





Figure 1. Profile of the radial electric and magnetic fields as funcions of p: the 
slope at the origin for / (see eqn. dT^l) ) obtained by our numerical shooting method 
is j3 ps 0.58319. The electric and magnetic fields in the plots have been rescaled 
appropriately. 



At this point, it is important to stress that the vortex solution for the present model 
does not exist in the absence of the Chern-Simons term. Indeed, a glimpse at equation 
02]) shows that the Chern-Simon coupling k induces a charge density proportional to the 
magnetic field b, which turns out to be crucial to source the vortex electric field. 



3.1.3. Energy, entropy and free energy: In the previous section we have shown the 
existence of vortex solution in the system. A natural question to ask is whether the 
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E v 



system will choose it as its ground state or not. In order to answer this question, we 
write the expression for the energy in the form 

where a is a UV cutoff, and L the size of the system. The vortex is regular at the 
origin, so we do not expect any singularity in the a — > limit. Taking into account the 
asymptotic behavior ( IT4l) one immediately finds that 

. , 7r|m| 4 / L\ , , 

E v « E core (R) + In ( -J , (16) 

where E core (R) is the finite contribution arising from the energy density integrated up 
to a distance R bigger than the core radius R Q . 

Since the vortex energy diverges logarithmically with the size of the system, the 
standard Kosterlitz-Thouless argument follows [H] : energetic considerations imply that 
at negligible temperatures the system will never choose the vortex configuration, but at 
finite temperatures, the choice of background is governed by the Helmholtz free energy 
F = E — T S and the entropic contribution could favor vortices for high enough T. In 
our vortex example, this entropic contribution arises from the number of possible places 
on the plane in which the vortex can sit, 

S = 2k B In (-) . (17) 

The similar logarithmic behavior for both the energy and entropy contributions will 
compete on the free energy, 

F = E - TS = 

= E core {R) + ^ In (|) - 2k B T\n (£ 

As a result, Kosterlitz and Thouless argued that a temperature should exists above 
which the system lowers its free energy by popping vortex-antivortex pairs out of the 
anisotropic vacuum [U],[T2]. This phenomenon is described in the literature as vortex 
unbinding and named topological phase transition. The critical temperature Tkt f° r 
this phase transition can be estimated from the F = condition, Taking into account 
the L dependence we find for the present case 

m 7r|m| 4 , . 

Below this critical temperature the system will be in a quasi-long-range ordered 
anisotropic phase. As soon as the system reaches Tkt it will be energetically favored to 
produce vortices and the quasi-long-range order will be destroyed. 



3.2. Hyperbolic plane 



Negative curvature spaces have been proposed as interesting setups to cure infrared 
divergences. The argument is pretty simple [8]: since the volume of space grows 
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Figure 2. Profile of the electric and magnetic fields in the hyperbolic plane as 
functions of p. For = 4, 2, 1.25 the corresponding slopes of the electric field at the 
origin, obtained by our numerical shooting method, are /3 « 2.01093,0.80291,0.30023 
respectively. 



exponentially with the distance to the origin, Gauss' law implies an exponential decay 
for massless fields. In this section we will analyze the modifications that arise on the 
vortex solution when we formulate the LCS model on the 2c?-hyperbolic plane[9]. 

As in the previous section, the vortex solution is easily found when writting the 
metric in polar coordinates, 

ds 2 = R 2 {dp 2 + sinh 2 p d0 2 ) , (20) 

here R is the curvature radius of the space. Inserting the static radial ansatz into 
the equations of motion ([5]) again with m 2 = — \m\ 2 leads to 

1 „, 1 



0. 



(21) 



tanhp" tanh z p' 

where f 2 = (R 2 + R 2 ) / 'R 2 . The existence of curvature in the 2-dimensional space results 
in an equation of motion now depending on a parameter f 2 which cannot be re-absorbed. 

The behavior for / near the origin coincides with that of flat space (j!3p . but the 
large distance behavior changes drastically. The relevant case for us is f 2 > 1 which 
results in a non-zero value at infinity for the electric field and an asymptotic radial 
profile given by 

2 



f(p) = V/o 2 - 1 1 



O 



(22) 



Notice that the exponential decay for the profile is independent of the parameters of the 
model. We have depicted in Fig.(j2]) the numerical solutions we find for different values 
of fl 

The vortex energy functional in the hyperbolic plane case takes the form 
n\m\*Rt r L . . , / „,o f 2 , n R 2 ff 



E 



{hyp) 



XR 2 



? 4 r L ( 
± J dpsmhplf' 2 + 



R„ 

_ _|_ 2—2- 

tanh p R 2 tanhp 



+ 



+ \f - %f + K 



(23) 



with a, L being respectively UV and IR cutoffs. The energy density behavior near 
the origin coincides with that of flat space and far away from the origin, although the 
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Figure 3. Profile of the electric and magnetic fields for the sphere as function of the 
azimuthal angle p (see eqn.([25jl). For = 3, 5, 10, the slope of the electric field at the 
origin are /3 « 1.91747,3.143601,6.090353. 



volume element grows exponentially fast at large distances, the electric field decays more 
rapidly (see ( 1221) ) and overwhelms the divergence. Taking the asymptotic behaviors into 
account we find that the vortex logarithmic divergence present in flat space is cured by 
the spatial negative curvature. We have therefore found that the n = 1 solution is a 
true soliton (finite energy) for the hyperbolic space case. The appropriate value for E Q 
results 

g = (-R 2 - R 2 - Rl)(R 2 - R 2 - 3R 2 Q ) . . 

2m*RtRi ' 

3.3. Sphere 

For completeness we will now consider the LCS model formulated on a sphere with 
metric 

ds 2 = R 2 [dp 2 + sin 2 p d6 2 ) , (25) 

here R is the radius of the sphere. Introducing the ansatz (j5J) into (jSJ), the equation of 
motion with m 2 = — \m\ 2 results in 

/" + r^—f - rV/ - / if ~ fo) = , (26) 
tan p tarn p v ' 

with f 2 = (R 2 — R^)/R 2 . To avoid singularities, we should demand the electric field to 
vanish at the north and south pole of the 2-sphere. The appropriate boundary conditions 
in the sphere case result 

e p ( P , 9)\ p=0 , n = => /(0) = /(tt) = . (27) 

The asymptotic behavior at the poles coincide with the flat space case and the 
configurations we obtain with these boundary conditions correspond to a vortex and 
an antivortex located at antipodes on the 2-sphere. This last fact is a consequence of 
(|2J): no net charge can be supported on a compact manifold. Therefore the flux lines 
emanating from the north pole should end on an equally opposite charge located in our 
case at the south pole. Some typical profiles are plotted in Fig. (j3J). 
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The energy for the vortex-antivortex configuration on the 2-sphere takes the form 
{sph) n\m\*Ri . / p 2R 2 ff> 1 4 

(28) 

In the sphere case the energy of the solution is finite due to the compactness of the 
spatial manifold and the regularity of the solution, there is no need for a compensating 
constant E . 



4. Closing Remarks 

We found vortex-like solutions to the Lifshitz-Chern-Simons theory. The presence of 
the Chern-Simons term is crucial to the existence of such solutions, since it sources the 
Gauss law (J2]). The logarithmically divergent energy suggests that a Kostelitz-Thouless 
transition may occur on the system. To fully clarify the nature of this phase transition 
a renormalization group analysis along the lines of [13J should be performed, this is 
currently underway and will be presented elsewhere. Unlike the U(l) global vortex 
case, extensions to higher winding are non trivial, due to the vector character of the 
electric field. 

Following a suggestion made in [8], we studied vortex solutions of the Lifshitz- 
Chern-Simons model formulated on the hyperbolic plane. We found, as expected, 
that the resulting configurations have finite energy. For completeness, we also explore 
Lifshitz-Chern-Simons vortex solutions on the sphere. In this last case, as in any 
compact manifold, the solution we found consisted in a vortex-antivortex pair. An 
open point is to study the stability of such solution. 

An important question is whether the flat space version of action jT]) at m — is 
the most general action describing the critical point. In principle, higher order terms 
like the introduced in [H] could be added, without breaking the z = 2 scaling symmetry. 
We will explore this issue in a separate publication. 
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Appendix A. Electromagnetic duality in Lifshitz systems 

In this appendix we will shortly review the electromagnetic duality mapping the k = 
action (CEJ) into az = 2 scalar field theory [I]. [17]. 

We start with the flat space action (CD) without the Chern-Simons deformation, 

S = J dtd 2 x L(d tai - diat) - i (^(V^) 2 + &)-\ (™ 2 (eif + -(e*) 4 )) • 

(A.l) 

The equations of motion following from the action are 

(k 2 V 2 - m 2 - Xe 2 )ei = d { a t - d t a h (A.2) 
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dte* = 0, (A.3) 
eijdjb + dta = 0. (A.4) 
Written in terms of the gauge invariant fields e« , b they read 

(ft 2 V 2 - m 2 )eijdiej - \e^di{e\ ej) = -d t b, (A. 5) 

die* = 0, (A.6) 
eijdjb + d t ei = 0. (A.7) 

The duality transformation is defined by solving Gauss law ( 1A.6I) as 

e* = e ij dj<f). (A.8) 

and then plugging (IA.80 into equation (1A.7I) . This implies that b can be written in terms 
of 4> as 

b = -d t <t>. (A.9) 

where a time dependent integration "constant" has been reabsorbed in <p. Replacing 
(fA~8l)-(fAT9l) into (fA~5l) leads to 

- (k 2 V 2 - m 2 )V 2 + A^((V0) 2 ^0) = «9 2 0. (A.10) 

This equation of motion for the scalar field can be derived from the action 

S[<j>] = J dtd 2 x l - (W) 2 - k 2 (V 2 0) 2 - m 2 (V0) 2 - \{V<P)^J (A.ll) 

The first two terms correspond to the z = 2 Lifshitz action for a scalar field. As it 
happens in its electromagnetic dual, the m 2 term drives the theory away from its z = 2 
fixed point into a z = 1 IR fixed point. Note that the gauge vector a i: a t decouples from 
and can be integrated out. 

In the presence of a CS term, the dualization described above leads to a non-local 
action for the scalar field [I]. An alternative non-local duality transformation leading 
to a local equation of motion for the scalar field can be found in [T8] . 



